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We obtain through a Matrix Product Ansatz the exact solution of the most general inhomogeneous 
spin chain with nearest neighbor interaction and with U(l) 2 and U(l) 3 symmetries. These models 
are related to the one loop mixing matrix of the Leigh-Strassler deformed N — 4 SYM theory, dual 
to type IIB string theory in the generalized Lunin-Maldacena backgrounds, in the sectors of two and 
three kinds of fields, respectively. The solutions presented here generalizes the results obtained by 
the author in a previous work for homogeneous spins chains with U(1) N symmetries in the sectors 
of N = 2 and N = 3. 
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INTRODUCTION 

String theory was first introduced in the 60's years, as a tentative to describe the big amount of mesons and hadrons 
discovered in particles accelerators. In this theory the particles are seen as different oscillations modes of the strings 
and it was able to describe some spectral characteristics of hadrons. Although this success, the original string theory 
is not able to explain many physical phenomenon mediated by strong interactions and nowadays the most successful 
theory to describe particle physics is the Quantum Cromodynamics (QCD). On the other hand, despite string theory 
has reemerged in the last three decades as a promising candidate for a quantum theory of all known interactions, its 
validity to describe nature has generated heated discussions since the theory is untestable due to the experimental 
impossibility to reach the tremendous energies found at the Plank scale. So, what is the actual relevance of string 
theory for real world physics? Either string theory is the correct theory to describe the world or not, remarkable 
resulties showed that string theory can be seen as a different formalism to quantum field theory. In QCD, it is not 
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possible to obtain a satisfactory quantitative description in the small energy regime, when the coupling constant is 
very strong. In this regime, numerical calculation on the lattice is the best tool for study physical models. On the 
other hand, it was observed by 't Hooft [T] that the theory is simplified when the number N of colors is very high, 
't Hooft derived a relationship between the topological structure of a Feynman graph and its N dependence. When 
1/N is interpreted as a coupling constant, an expansion in 1/N is similar to an expansion in a generic interacting 
string theory, resulting in a relation between strings and planar diagrams. If the N = 3 can be regarded as a large N, 
it explain why the string models of 60 's years was able to give the right relation for some spectral characteristics of 
hadrons. More recently a remarkable result arouses new interest in the duality between string theory and quantum 
field theory. Maldacena conjectured that IIB string theory on the curved background AdS^ x S 5 (anti-de Sitter and 
sphere spaces) should be equivalent to N — 4 Super Yang-Mills (SYM) [3j|4j|5]. This conjecture, AdS/CFT, relates 
operators, states, correlation functions and dynamics of both theories. One of the most important results of this 
conjecture predicts that the spectrum of scaling dimension operator of gauge invariant operators, in the conformal 
field theory, should coincide with the spectrum of energies E of string states. Furthermore, this correspondence relates 
the weak coupling constant regime, in the gauge theory, with the strong coupling constant ones, in the string theory. 

The Maldacena' s conjecture need yet to be fully proved and since the discovery of the relation between the 
planar dilatation operator of the Af = 4 SYM with an integrable so(6) quantum spin chains [2], integrability has 
played a prominent role in the exploration of the Maldacena's correspondence. The study of the planar dilatation 
operator's integrability is very important because it not only enable us to test the Maldacena's correspondence as 
it is an generator of nontrivial integrable models. Exactly solvable models are of interest in high energy physics, 
condensed matter physics, statistical mechanics and mathematics since the pioneering work of Hans Bethe [5] (see, 
e.g., |7j |9j [10] for reviews). According to this ansatz the amplitudes of the eigenfunction are expressed by a 
nonlinear combination of properly defined plane waves. On the other hand, in the last two decades several different 
ansatz were introduced in the literature under the general name of matrix product ansatz (MPA) . The first formulation 
was done for the description of the ground-state eigenfunction of some special non-integrable quantum chains, the so 
called valence-bond solid models [TTJ [T5J [T21 Q3] . The MPA becomes also a successful tool for the exact calculation 
of the stationary probability distribution of some stochastic one dimensional systems [15] [THJ [17] . An extension of 
this last MPA, called dynamical MPA was introduced in [T5J [H5] and extended in [20 . This last ansatz gives the 
time-dependent probability distribution for some exact integrable systems. The MPA we are going to use in this 
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paper was introduced in [5TJ [53] [22]. This ansatz was applied with success in the evaluation of the spectra of 
several integrable quantum Hamiltonians [2TJ [22j [23] , transfer matrices [25] |26j [27] and the time-evolution operator of 
stochastic systems 24J. According to this ansatz, the amplitudes of the eigenfunctions are given in terms of a product 
of matrices where the matrices obey appropriated algebraic relations. In the case of the Bethe ansatz the spectral 
parameters and the amplitudes of the plane waves are fixed, apart from a normalization constant, by the eigenvalue 
equation of the Hamiltonian or transfer matrix. On the other hand, in the MPA the eigenvalue equation fixes the 
commutation relations of the matrices defining the ansatz. In such case the spectrum of the Hamiltonian or transfer 
matrix, and the corresponding eigenfunctions, can be computed in a purely algebraic way. 

In the present paper we obtain through a MPA the exact solution of the most general inhomogeneous spin chain with 
nearest neighbor interaction and with U(l) 2 and U(l) 3 symmetries. This model is related to the one loop dilatation 
operator in deformed Lunin-Maldacena backgrounds |28] and conformal field theories with deformations [29J. The 
solutions presented here generalizes the results obtained for homogeneous spins chains with U(1) N symmetries in the 
sectors of N = 2 and N = 3 [30] [31]. In this model the coupling interaction between neighbor sites are not a constant 
as it is in the homogeneous model studied in [3U1 13T] , 

THE INHOMOGENEOUS SPIN CHAINS 

We consider here the most general inhomogeneous 3-state spin chain with nearest neighbor interaction, periodic 
boundary condition and U(l) 3 symmetry. The C/(l) 3 symmetry imply that the Hamiltonian describing the time 
evolution of this spin chain conserves the number of states of each type. By denoting the basis of states at a given 
site as \a) (a = 1, 2, 3), the Hamiltonian in a periodic lattice with L sites takes the form 

H = t,l E r^ a (jj + i)E^Ef + {+ £ r^(j,i + i)^^), (i) 

2=1 \a#3=l a, 0=1 J 

where E a @ are 3x3 Wyel matrices with elements (i? Q/3 )._. — Sa^Spj (a,f3 — 1,2,3). While the first term in 
the right hand side of ([I]) acts over neighbor sites exchanging its configuration \a)j <8> — ► ® with 
rate T^^(j,j + 1), the second one is a diagonal operator with weig ht T^(j,j + 1). The Hamiltonian |j} m a 
particular case (T"g(J,j + 1) = T"^ constant) contain the homogeneous integrable spin-1 chain, related to the one 
loop dilatation operator in deformed Lunin-Maldacena backgrounds, studied by Freyhult, Kristjansen and Mansson 



[3"T] and generalized in [3D]. The eigenfunctions for (|T|) can be construct as 

|*t>= ]T ^ ai - aL \ai,-,a L ) (ay = 1,2,3), (2) 

a lt ...,a L 

where the symbol (*) in the sum denotes the restriction to the sets {a\, with the same number n a of spins in 

configuration a and \]/ Q i - Q £ is a constant. 

THE MPA 

In order to formulate a MPA for the Hamiltonian ([TJ , we make a one-to-one correspondence between configurations 
of spins and products of abstract matrices. This matrix product is construct by making a correspondence between a 
site j in the chain with spin configuration a — 1, 2, 3 and a matrix Aj . Our MPA asserts that the components of the 
amplitude of the eigenfunction \J/ a i •••<*£ i n ^ are obtained by associating them to a products of these matrices Aj** . 
Actually are abstract operators with an associative product. A well defined eigenfunction is obtained, apart 
from a normalization factor, if all the amplitudes are related uniquely. Equivalently, in the subset of words (products 
of matrices) in the algebra containing n a {n\ + 112 + ^13 — L) matrices Aj there exists only a single independent 
word. The relation between any two words gives the ratio between the corresponding amplitudes of the components 
of the eigenfunction To formulated the ansatz we can choose any uniform operation on the matrix products 

that gives a non-zero scalar to make a map between the amplitudes \[/ Q i - Q i i n ([2| and a matrix product [55] : 

V ai - aL ^ A[ ai) A { 2 a2) ■■■A<£ l) (ay = 1,2,3). (3) 

It is obvious that the 3 states \a...a) (a = 1,2,3) are all eigenstates of the Hamiltonian ([!]). In the following we 
shall choose |l...l) as our reference state. The Hamiltonian ([!]) when applied to the components of the eigenfunction 
([2]) where we do not have spins configurations |a) (a = 2, 3) at nearest neighbor sites and at boundary positions give 
us the constraints, for the amplitudes yr a i— a L ([3j) ; 

EnA 11 " 1 ^ 1 '! 12-11-1 • • • A x aj -^A Xj ~ X3 - 1 ~ 1 A < - x a ^ A Xj + l ~ x o- 1 A ( ^ lj+1 ' ) ■ ■ ■ A x an) A L ~ Xn 

n 

= [r^.C?, j + ^A^A^A^ 1 - 1 • ..^Hl^r^-i^^^^+i-^^jti) . . . A {*n) A L- Xn 
3=1 

+ ai \ x (j,j + i^a^^a^a**-^- 1 ■ • . i^-o^i-^-i^il i^+i-^-^^+o . .. a m a l-x„ 

L 

+ E [ T Z 1(M + 1) + rj Zl 1 + !)] A'^Ai^A"-**- 1 ■ ■ ■ A X ^A L ~ X ~ (a, = 1, 2, 3), (4) 
1=1 



where e n is the energy of the eigenfunction (j^J), A = vli 1 '', n = ri2 + ^3, and x%, -;X n are the position in the spin chain 
where we have a state configuration \a ^ 1). A convenient solution of this last equation is obtained by identifying 
the matrices A x {ot — 2,3) as composed by spectral-parameter-dependent matrices. The distinguibility of states 
configurations allows two types of solutions. The standard solution is obtained if each of the matrices A^ (a = 2, 3) 
is composed of n = n 2 + 113 spectral parameter dependent matrices 22 , 24, 30]. A second class of solutions is obtained 
if matrices Ax (a — 2, 3) with different a value are composed of by distinct sets of spectral parameters matrices 
[22]. [30] . Here we will consider only the standard solution but our model can be easily extended to the second class 
problem. In the present case, the matrices A x (a — 2, 3) can be written in terms of the matrix A and n = n 2 + 113 
spectral parameter dependent matrices A " k [3"9"| : 

n 

4 Q) =E43A (« = 2,3), (5) 
i=i 

where the matrices A^fj, . satisfy the following commutation relations with the matrix A: 

A { S j A = g a (x,x + l)e^AA x ^ hk ., (a — 2, N), (j = (6) 

the parameters kj (J = 1, ...,n) are in general complex numbers unknown a priori, and g a (x,x + 1) is a constant. 
The energy e„ is obtained by inserting ^ in Q, by using (JsJ) and imposing that e n is a symmetric function on the 
spectral parameters [30] 

n 3 

e n = £ (r? + T\ \e- ik ') + ]T n a (T\° + T« \) + (L- 2n)T{ {, (7) 



where we need to impose 



rjj = r^(x,.x + 1) 
r«i(x,x + i) m 



? a (x,x+l) = ^ aU Jl ■ ^ = ; (a = 2,3) (8) 



and 



rifta:,a! + l)=rii, r} "(a;, z + 1) = T\ «, T% {(x, x + 1) = T% { (a = 2,3), (9) 

where I 1 ^, r£ ?, rj g and } are constants. 

The relations coming from the eigenvalue equation for configurations where we have two spins configurations \a) 



(a = 2, N) at nearest neighbor sites and are not located at boundary positions given us 



n 

E [ r 2 1 + r ? WV»+ k * + (T a a \ + r{° - T\ \ - TZZ(x,x + l))e^ 



i,i=2 



E [ri? + r? £e^+*«> + (rs 1 + rj g - rj J - r» a: + i))e** 



\p riq( x,X+ l) T ,^a / 



Apt) _ 
iki AP) A a ) 



(10) 



j,l=2 1 1 M ' x + 

where we have used ([3|, ([5|-([9| and y — 1, ...,L. The relations (10 1 should be satisfied for all x,y — 1, ...L. This is 



possible if we have constants ™, g, ^ and if 



rfl(x,x + i) 



r« «(x, x + 1) = r« gfo * + 1) = v a a % r£ * + 1) = 



./3 a. 



Finally, the relations jlo| fix the algebraic relations among the matrices {a = 2, 3): 



(ii) 



(12) 



a',/3'=2 



Relations (|6j) and (12 1 define completely the algebra whose structural constants are the S*-matrix of the spin-| model 
[3TJ [52 ED] . Since the several components of the wavefunction should be uniquely related, the above algebra should 
be associative. This associativity implies that the above ^-matrix should satisfy the Yang-Baxter relations [71 132). 
which is indeed the case |30j . The components of the wavefunction corresponding to the configurations where we have 
three or four particles in next-neigbouring sites would give in principle new relations involving three or four matrices 
■^xk - These new relations are however consequences of the above relations (6 1 and ( j 1 2 j ) - It is important to mention 
that in the sector ri2 = or 713 = (symmetry U{\) 2 ) the Hamiltonian |l]) reduces to a inhomogeneous version of the 
well known asymmetric XXZ model |33j . 

In order to complete our solutions through the MPA Q we should fix the spectral parameters, or momenta, 
ki,...,k n . Theses parameter are fixed from the configurations where we have a spin configuration a = 2, 3 at 
boundary positions (x\ = 1 or x n = L). By using the algebraic relations ([5]), (|6| and (12 1 we obtain the relation 



N 



E (ai,---,a n 

a' x ,...,a' n =2 



(13) 



where we have used the identity (see [30] ) 



3 J 



(14) 



and 



<<*!,..., a n \T^\a[,...,a' n )= £ {s£ $ (fc x , A*) • • • S% ^ k s ) ■ ■ ■ S% $ (k n , (15) 

a",...,a^ 

where 4>{ot") = X\ x —\ 9a'{( x , x + 1), is a (2)" x (2) "-dimensional transfer matrix of an inhomogeneous vertex model 



(inhomogcncities {k{\) with Boltzmann weights given by the S'-matrix elements defined in ( 12 1. The model is defined 



on a cylinder of perimeter n with a seam along its axis producing the twisted boundary condition 

^v5 +1 (^) = Syikkn.k^a'i). (16) 



Finally relation (13) with (16 1 give us the constraints for the spectral parameters: 

e -ik j L =A (n) {k , ){h}) (i=1) ... jn)) (17) 



where A( n )(fe i , {&;}) are the eigenvalues of the transfer matrix (U5|. The condition (17k leads to the problem of 
evaluation the eigenvalues of the inhomogeneous transfer matrix ( |15[ ) . This can be done through the algebraic Bethe 
ansatz [34 or the coordinate Bethe ansatz (see [35] and [35] for example). 

DISCUSSION AND CONCLUSIONS 

We solve through a MPA the most general inhomogeneous 3-state spin chain with U(l) 3 symmetry and nearest 
neighbor interaction. We found that the coupling constants in ([lj should satisfy the constraints ([8]), Q and (111 in 



order to make the Hamiltonian (JT|) integrable. It is important to mention that different from the homogeneous spin 
chain, where the eigenstates Q are also eigenstates of the translation operator due to the periodic boundary condition, 
in the inhomogeneous model the eigenstates do not have a defined momentum. It happens due to the inhomogeneities 
that broke the translational invariance of the system. The study of this new model can be of interest in the context 
of both AdS/CFT and condensed matter physics since it is related to the one loop dilatation operator in deformed 
Lunin-Maldacena backgrounds [28 , conformal field theories with deformations [29J and inhomogeneous spin chains. 
Another quite interesting problem for the future concerns the formulation of the MPA for the case where we have 
open boundary conditions, as well as for quantum chains with no global conservation law such as the XYZ model, the 
8- vertex model or the case where the quantum chains are defined on open lattices with non-diagonal boundary fields. 
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